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Abstract 

The purpose of the present work is to apply the method recently developed in reference to 
the spin-1 Ising chain, showing how to obtain analytical /3-expansions of thermodynamical func- 
tions through this formalism. In this method, we do not solve any transfer matrix-like equations. 
A comparison between the /3-expansions of the specific heat and the magnetic susceptibility for 
the s = ^ and s = 1 one-dimensional Ising models is presented. We show that those expansions 
have poorer convergence when the auxiliary function of the model has singularities. 



1 Introduction 



One dimensional (ID) lattice models have been widely studied for decades, motivated, for 
example, by the possibility of obtaining physical insights about more realistic theories and of their 
utilization as good laboratories for testing the applicability of new methods. 

Several instances of chain models have been intensively investigated, mainly by virtue of their 
relevance in connection with Condensed Matter phenomena; among them, we can mention the one- 
dimensional (ID) Ising, the spin-1/2 XXZ Heisenberg^] and the Hubbard^, Q models. Actually, the 
chain models mentioned above belong to the class of exactly solvable models, due to their integrability 
property 0. Thus, the algebraic Bethe ansafe method has been widely applied in order to determine 
their thermodynamical functions by means of non-linear integral equations. 

Recently Rojas et al. showed how to obtain a closed analytical expansion for the grand 
potential in the thermodynamic limit for any translationally invariant chain model from the high 
temperature expansion of the cumulant method^. Another point that we can stress is that the 
method developed in ref. can also be extended to models including nearest interacting neighbors 
(e.g. frustrated quantum Heisenberg models with spin-Sp, 0). For such quantum systems, the 
calculation of higher order terms in the /3-expansion of thermodynamical functions allows one to 
study their properties at lower temperatures 0. Namely, very good results have been obtained 
applying the method to the (integrable) spin-1/2 XXZ chain model[0]. Actually, we should stress 
that the applicability of such method extends further to either integrable or non-integrable models 
(the XXZ with spin-l-case being an example of the latter case), provided that the hamiltonian 
presents translational invariance, nearest-neighbor interaction, and periodic boundary conditions. 

In this paper we intend to obtain an analytical /?-expansion of the Helmholtz free energy for the 
spin-1 Ising model through the application of the results of reference Using the fact that this 
model is exactly solvable, the correctness of such expansion can be checked. The derivation of the 
Helmholtz free energy for any unidimensional chain model in the thermodynamical limit is obtained 
from the auxiliary function v'(A), the latter representing the actual mathematical improvement 
in this approach, in relation to the literature (cf.^]). Moreover, the present paper is a first step 
towards the application of such approach to ID Ising models with arbitrary spin|§], as well as to 
more interesting models, such as the spin-1 XXZ Heisenberg chain 

The importance and usefulness of the /3-expansion of thermodynamic functions relies on the 
perspective of obtaining information for lower temperatures, as terms of higher order in the expansion 
are taken into account. Hence, knowing how fast or how slow such expansion converges is crucial to its 
applicability. We want to clarify how the characteristics (namely, the singularities) of the auxiliary 
function ^{X) impact on the convergence of the /3-expansion of any thermodynamic function. In 
order to do so, we have studied the /3-expansion of the thermodynamic functions for the spin-1 and 
spin-1/2 Ising models. At this point, we should stress once more that our /5-expansion results are 
analytic, non-perturbative and valid, in principle, for arbitrary values of the parameters. 

In section 2 we present a summary of the results derived in reference ||l|. In section 3 we obtain 
an analytical expression for the auxiliary function ^{X) for the spin-1 Ising model for finite values of 
13] in section 4 we discuss the rapidity of the convergence of the /3-expansion of two thermodynamical 
functions, namely, the specific heat and the correlation between the z-component of spin between 
nearest neighbors, when the auxiliary function ^p{\) is/is not singular at A = 1. For the spin-1 Ising 
model, we restrict ourselves to the case /i = 0. In section 5 we present our conclusions. 



2 A Survey of the Calculation of the Grand Potential of a 
Chain Model 



In reference ^ Rojas et al. obtained, in the thermodynamic hmit, a closed expression for the 
grand potential of any translationally invariant chain model with periodic boundary condition from 
an auxiliary function. Here we present a survey of the method and we refer to reference [||] for further 
details. 

Let us consider a one-dimensional regular lattice (a periodic chain) with sites, so that the 

N 

Hilbert space of the chain model is simply THS ' = ®T-i, Ti being the irreducible representation at one 
site, including all its degrees to freedom. The dimension of this Hilbert space is dimTi*^^^ = tr7v(]l). 
The notation tr^r means the trace over all sites and their internal degrees of freedom, e.g. spin. 
The grand canonical partition function of a quantum system in the chain with sites is given 

by 



ZAr(/?,/i) = tr^(e-^^), (1) 

where K = EI — ^uN, with /i being the chemical potential and N being an operator that commutes 
with the Hamiltonian of the system. 

Let A be any operator that acts on 7i*^^^^ where M ^ N. We define (A) = *4 m ffl^ ' ■^^^ ^^-^ 
dimension of TC^^\ We call (A) the normalized trace of operator A. 

Using the definition of normalized trace, eq.([l|) becomes the expansion of jj) around /? = 0, 



= tr;v(]l){l + 



n=l 

In reference we showed that the coefficients (K") in eq.(|I]) can be written, for any translationally 
invariant hamiltonian with interactions between first neighbors and for arbitrary n as 

r=l m=r ^ ^ 

The notation [n, A^] means the min(n, A^) and is defined by 

n m r 

KS-E"E"nKS' (4) 

{ui} {mi} 3=1 

n 



means the restriction: ^rij = n and 7^ for i = 1, 2, .., m. We also use the notation 

{rii} 

{ui} = {rii, 712, ■ ■ ■ , and {nii} = {mi, m2, . . . , mm}- The function K^"2i is defined as 



in) 



{■Hi} 1=1 



(5) 



and each term on the r.h.s. of eq.@ corresponds to the g-trace of an open connected sub-chain. 

(n) 

In the definition of the function K\ ^ we have the g-trace which means 

\ n,\n2\ ...nj. / ' n\ l^n,n+i ' ^.2.2+1 ' • • • ' ^w.w+iJ^ l^^j 

where ^^j^ rii = n with rij 7^ and the indices ik, k = l..m are distinct among themselves. The nota- 
tion (P(Kj^^j^+i , Kj2,i2+i , • • • , ^im,hn+i)) fspresents all the distinct permutations of the m operators 

We show in reference that the grand potential per site VV(/5,/x), in the thermodynamic limit 
is written as 

/x) = {ln(tri(l)) + ln(l + 0} • (7) 

where 



d" ( Axr+' \ 



A=l 



with A being a parameter, the auxiliary function ip{X) is equal to 



m,=l 



(9a) 



and 



7i=m 

From eqs.(|7|)-(|9li|) we see that the grand potential per site, in the thermodynamic limit, can be 
derived only from the open connected sub-chains. The weight of each sub-chain in the /3-expansion 
of W(/3, /i) is already presented in eq.(|^. 

3 The Helmholtz Free Energy of the Spin-1 Ising model 

The hamiltonian of the spin-1 Ising model (with single-ion anisotropy) is: 

N 

M = J2 ^S;S]^^ - hS^ + D{S'^)\ (10) 
i=i 

where h is the external magnetic field in the z-direction, while the parameters A and D are the 
exchange and single-ion anisotropy, respectively. The chain has sites, and periodic boundary 
conditions are assumed. 



It is worth noticing that the partition function for the model above may be obtained from the 
transfer matrix approach by calculating the eigenvalues of the hamiltonian (p^. These eigenvalues 
are solutions of a polynomial of third degree. Nevertheless, for arbitrary values of the parameters A, 
h and D, such an approach only gives numerical solutions for the thermodynamic functions, which 
allow us to verify the applicability of the method presented in reference fTj] to one exactly solvable 
limiting case of a non-integrable model 0]. The validity of the present approach has already been 
verified for the integrable XXZ Heisenberg spin-1/2 modeljjl]] and its limiting cases |1I0[ . 



In order to obtain the Helmholtz free energy, we must calculate at first its functions H^"^. Due 
to the fact that all terms that contribute to the hamiltonian of the Ising model commute among 
themselves, we may write 

" m Tjni 



{m} 1=1 



<-=E(n^)' with m<n. (11^ 



The sums in Tm (see eq.(^)), including the restricted sums over the indices {rrii}, are more tractable 
if we recognize that they can be substituted, in the thermodynamic limit, by m independent sums 
with each index {rrii} varying from 1 to oo[|l|. For this model, becomes 

ni,n2,--- ,nm=l i=l 

where = AS'fS'f^.^ + i?f , with i?f = —hS^ + -0(5"/)^. The two operators in commute; 

consequently, we can apply Newton's multinomial formula to obtain the expansion of H'^l^i, 

^+1 = E f?) (A)"'-^^(5?)"-^^(5n^15^i)"'-^% (13) 

where multinomial coefficients. After averaging over the space sites from i = 1 up to 

i = m + 1, we obtain 

i=l \ n^=l j,=0 ^ 

x[l + (-l)'^™-^-+5,„,nJ. (14) 

It is tedious but simple to show that 



oo 

Um^. ^ \J^J 

"■m — J- Jm — 'J 

x[l + + s.^^r^J = (-l)"™--^--ai + h, (15) 



1 V i-P) 
3 



where the constants ai and bi have been defined as 



ai = p'^ + q +r and bi = p + + (16) 

and 



Performing the sums in the product on the r.h.s. of eq. ([l4D , we get a recursive solution 



-L r 



m 



where 



ai = hi-ip^ + ai^iq and bi = bi^ip +ai^iq~^, = 2, 3, ■ ■ ■ , m. (19) 

Eq.(^) gives us the relation between the auxihary function (p{X) and which, for this particular 
model, is 



m=l »n=l 

Defining (p{X) = 0a(A) + 0b(A) with 

(21) 

m=l m=l 

we may use the relations (|T9t) to rewrite 0a(A) and ^^(A) as 



and 



0a(A) = II + ^0,(A) + ^0„(A) (22a) 



06(A) = ^ + ^0,(A) + ^0„(A) (22b) 



which yields 



1 [p^ + p + q'^ + q + r+ + r )A + 2(p+g+— p q )+r^(p^ — q )+r (g^— p ) 

'^(^'' = ^ — { ~ , ~\\ — nrr — • (23) 

3 A^ — (p+gjA — [p+q+ — p q ) 

The Helmholtz free energy of the Ising model is then obtained by substituting eq. (|23|) in eq.(0), that 
is. 



w(.).-i{M3)..(:.|:^(^) 

and from it we derive a /3-expansion of with arbitrary value of n through an algebraic compu- 

tational language such as MAPLE. 

Depending on the values of A, h and D, the function v'(A) may be singular at A = 1 for some 
values of (3. This is also true for powers of this function at A = 1 which contribute to eg. (p^ . 
However, we know that there is no phase transition at finite j3 in any one- dimensional Ising model. 
Consequently, those singularities in ip{l) (i.e., ip{\ = 1)) must be non-physical. We point out that a 
similar behavior is shown by the auxiliary function for the spin-1/2 Ising model (see reference [p!0[]). 
A question yet to be answered is whether the presence of non-physical singularities in v'(l) could 
indicate a poorer convergence of high temperature expansions of the thermodynamical quantities. 
In the next section we compare the high temperature expansions for the Ising models with spin-1/2 
and spin-1, when their respective auxiliar functions have a singularity at A = 1, and when they do 
not. 




4 Comparison of the High Temperature Expansions of 
s = I and s = 1 Ising models 

The auxiliary function v^(A) appears for the first time in the literature in our reference 
Therefore calculating it for exactly solvable models is useful to the understanding of its properties, 
including the impact of its singularities on the convergence of the /3-expansion of thermodynamic 
quantitities. 

The function (pi{X) of the spin-1 Ising model (see eq.(^)) can be written as 



^i(A) = ^^ + 3^, (25a) 
where X± are the roots of the polynomial of second degree in A in the denominator of eq. (|23|) , namely. 



Ah 



+ p ± a/ {q- - p-)2 + Aq- 



■p-\ 



(25b) 



The constants A{f5) and -B(/3) are easily obtained substituting eq. ( |25a| ) in eq. 
For the case h = 0, eqs. (pSD reduce to 



v^r^(A) 



A{(3) 



A - (p+ +p^)' 



(26a) 



with 



(26b) 



For h = 0, the type of A-dependence in eq.( p6aD gives rise to an alternating series to the Helmholtz 
free enerev, in the variable A(B)/(1 — X,V. 



For —D < A < D and D > 0, the function ip^^\l) has no singularity; otherwise, there is one real 
and positive value of (5 for which it is singular. 

The auxiliary function ( |26a| ) has the same A-dependence as the auxiliary function of the spin-1/2 
Ising model, that was summed up in reference |]T0|- Substituting eg. ( |26a| ) in eq.(^), and following 
the steps described in |jlO|, we obtain the Helmholtz free energy of the spin-1 Ising model at = 0, 



>Vi(/3) = -iln 



(27) 



valid for finite values of (3. This result coincides with the one derived from the transfer matrix 
approach, as well as with those obtained by numerical analysis. 



We derived in reference |]T0| the auxiliary function v9i(A) for the spin-1/2 Ising model. For 
arbitrary value of the external magnetic field h we got 



g2/3h I g2/3(h+A) _ 2 

nW = 1 3^^^ I- (2^ 

Only for h ^ Q his function has one real and positive value of (3 where v^i(l) is singular. The 
Helmholtz free energy of the spin-1/2 Ising model, for arbitrary value of h is||ll|| 



(29) 

that is also valid for finite values of j3. 

In both equations ( pOj) and ( p7| ) , the respective Helmholtz free energies are /3-expansions of infinite 
range. We are interested in the impact of the singularities of (^(A) on the rapidity of convergence 
of the (3- expansion of thermodynamical quantities. For both spin=l and spin-1/2 models, we have 
taken the specific heat and the z-component of spin correlation < S^Sfj^^ > between nearest 
neighbors as examples, and have chosen two suitables sets of values for the parameters (A, D, h) so 
that the singular and non-singular cases can be compared, as far as the rapidity of convergence 
of /5-expansions is concerned. In what follows, approximate curves will refer to those obtained by 
truncation of the /5-expansion at 80th order in /3, in contrast to the exact curves. 

For the spin-1 Ising model, we have taken the first set of conditions to be (A = 1,D = l,h = 0), 
which yields a non-singular ^^[^■'(l); by promoting a slight variation in D, a second set {A = 1, D = 
0.99, h = 0) can be defined, yielding a singular (ff'\l). For both sets, exact and approximate curves 
of the specific heat Cy are shown in figure 1 and 2. We expect that both exact curves be slightly apart 
from each other, in some finite interval of sufficiently small /3 > 0. They differ, indeed, by a relative 
error inferior to 1%, for < f3 < 2.5. If we compare the behavior of the corresponding approximate 
curves, however, we observe that convergence in the singular ip case is much worse than that of the 
non-singular case. In the non-singular case, the approximate and exact curves differ by less than 1% 
for < (3 < 1.16, whereas in the singular case they differ by less than 1% for a much smaller interval 
< /3 < 0.54. Similar behavior can be observed in the correlation function < S^Sf_^_i > (see figures 
3 and 4). Exact curves in both singular and non-singular cases differ by less than 2% in the interval 
< j3 < 1.5; however, in the non-singular case the approximate and exact curves differ by less than 
1% in < /5 < 1.23, whereas in the singular case this interval is much shorter, namely, < (3 < 0.58. 



>V.(/3) = --ln 



/3A 

' 4 cosh(/9/i) + 



e '^2 sinh^(/3/i) -|- 



e 2 



Similar discussion can be carried out for the spin-1/2 Ising model. Here, the conditions (A = 
— 1, D = 1, h = 0) and (A = —1, D = 1, h = 0.01) yield non-singular and singular cases, respectively. 
(Observe the slight variation in h, only.) For both sets, exact and approximate curves of the specific 
heat C„ are shown in figure 5 and 6. The exact curves differ by less than 1% in the interval < 
f5 < 2.6. In the non-singular case, the approximate and exact curves differ by less than 1% for 
< f3 < 5.2, whereas in the singular case the same difference holds for a much smaller interval 
< /3 < 1.8. The same behavior is also exhibited by the correlation function < SfS^j^^ > (see figures 
7 and 8). Exact curves in both singular and non-singular cases differ by less than 0.05% in the 
interval < f3 < 1.8; however, in the non-singular case the approximate and exact curves differ by 
less than 1% in < /? < 1.23, whereas in the singular case this interval is much shorter, < /3 < 0.58. 

Those results indicate that the presence of singularities in the auxiliary function ip yields poor 
convergence of /3-expansions of thermodynamical functions. 



5 Conclusions 

In reference [0 we obtained, in the thermodynamical limit, a closed analytical series for the 
grand potential for any unidimensional chain model with periodic boundary conditions, from the 
high-temperature expansion of the cumulant method In order to derive the grand potential, 
we must calculate the auxiliary function v^(A) for the particular model of interest. Recently, the 
method presented by Rojas et al. was applied successfully to the spin-1/2 XXZ Heinsenberg model 
|l[] and its limiting cases|10|. In the present work we consider the exactly solvable spin-1 Ising model. 



whose numerical solutions can be easily attained and used in the verification of the correctness of the 
results derived from reference |^, as well as in the study of the properties of the auxiliary function 
(/9(A) of this model. From eqs. ( p3|) and (plj) we obtain the /^-expansion for the Helmholtz free 
energy VV(/3) for arbitrary values of the parameters A, D and h. The analyticity of this /9-expansion 
allows its use as input to the perturbative study of thermodynamical properties of uniaxial Ising- 
like models. For one set of values of the parameters, and n = 40 (where n is the leading order of 
the /3-expansion) we plotted those expansions and confirmed that our results match the numerical 
results for this thermodynamic function, in an interval of /5, including situations in which h 0. For 
h = we are able to sum the terms in the expansion ( |2^ of VVi(/9) for arbitrary values of A and 
D at any finite (3. In this case, the auxiliary function ipi\X) is given by eq. (^), and it has no 
singularity at A = 1 for —D < A < D and D > 0; otherwise, there is a value of (3 for which (p^i\X) 
is singular. Such singularity for the spin-1 model is of the same type as that of spin-1/2 Ising model 
||10|| . In both models the singularities are non-physical, since no unidimensional Ising model has 
phase transitions at finite (3. In order to see if the existence of a singularity in the auxiliary functions 
of both models could influence the rapidity of convergence of the /3-expansion of the thermodynamic 
functions, for each model we considered two distinct sets of parameter values: in one of them the 
respective auxiliary function is singular, and in the other one it has no singularities, at A = 1. For 
each model, the two sets differ by the value of only one parameter, and the difference in value is very 
small, in such a way that the results for those sets could be related to one another by perturbation 
theory, at least in a finite range of (3. 

In the /5-expansion of thermodynamic functions we kept n = 80 for both models and both 
cases (singular and non-singular (p's , spin-1/2 and spin-1). These results strongly indicate that the 
presence of singularities in the auxiliary function allows one to infer that the /3-expansions will have 
poor convergence, even though the series ( P3[ ) has infinite range in (3. 

Our analytical results at h = can be used as the starting point for analytical perturbative 
expansions for h ^ 0, and also for other models. 
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Figure 1: The full hue is the exact result of 
the specific heat for A = 1 and D — 1. The 
dotted line is its expansion up to order n=80 
in /? e [0,1.22]. 
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Figure 3: The full hne is the exact result 
of the correlation between nearest neighbors 
for A = 1 and D = 1. The dotted line is its 
expansion up to order n=80 in (3 e [0, 1.3]. 
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Figure 2: The full line is the exact result of 
the specific heat for A = 1 and D = 0.99. 
The dotted line is its expansion up to order 
n=80 in /3 e [0,0.57]. 



Correlation between nearest 
neighbors (s = 1) 

0.2 0.4 0.6 




Figure 4: The full line is the exact result of 
the correlation between nearest neighbors for 
A = 1 and i:) = 0.99. The dotted line is its 
expansion up to order n=80 in /3 e [0, 0.62]. 
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Figure 5: The full line is the exact result of 
the specific heat for A = — 1 and D = 1. The 
dotted line is its expansion up to order n=80 
in /5 e [0,5.85]. 
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Figure 7: The full line is the exact result of 
the correlation between nearest neighbors for 
A = — 1 and D = 1. The dotted line is its 
expansion up to order n=80 in (3 e [0, 6.3]. 
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Figure 6: The full hne is the exact result of 
the specific heat for A = —1 and D — 1, but 
for h — 0.01. The dotted line is its expansion 
up to order n=80 in /? e [0, 1.85]. 
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Figure 8: The full line is the exact result of 
the correlation between nearest neighbors for 

A = -1 and D = 1 but for h = 0.01. The 
dotted line is its expansion up to order n=80 
in /5 e [0,1.9]. 



